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Abstract. The availaible data on the D° — » K^TC + n" decays are analysed in the frame- 
work of the quasi two-body QCD factorization approximation. The annihilation, via W- 
exchange, amplitudes are added to the weak-decay tree amplitudes. The doubly Cabibbo 
suppressed parts of the amplitudes are also considered. The strong interactions between 
the kaon-pion and pion-pion pairs in the S -, P- and £)-nnal states are described in terms 
of the corresponding form factors. The kaon-pion or pion-pion scalar and vector form fac- 
tors are constrained by other experimental data. Unitarity, analyticity and chiral symmetry 
are also used to obtain their functional forms. We go through a minimization procedure 
to reproduce the K°n~, K°n + and 7t + n~ effective mass projections of the Dalitz plot dis- 
f> | tributions. The large number (27) of non-zero amplitudes leads to a large number of pa- 

rameters. The resulting model distributions and branching fractions are compared to the 
accurate Belle Collaboration data. 

1 Introduction 

Dalitz-plot time dependent amplitude analyses have been recently performed tl 1 121 for the CP self 
conjugate D° meson decays into K^n + 7T~. These studies have allowed a direct measure of the D°-D° 
mixing parameters, the knowledge of which could show the presence of new physics contribution 
beyond the standard model. Studies of — > D^K*, D — > K^tt + t:^ decays, in which the interference 
between D° and D° mesons was used to measure the Cabibbo-Kobayashi-Maskawa angle y, have also 
been accomplished 03141 . A good understanding of the final state interactions between mesons in the 
D° decays into K^n + n" is essential in order to reduce errors of the D° - D° mixing parameters and 
the measured values of the angle y. A very rich resonance spectrum seen in the Dalitz plots is a direct 
signal of the complexity of the strong meson interactions. Using these high statistics data theoretical 
models of the decay amplitudes can be tested. 

The experimental analyses Ml 121 relied mainly on application of an isobar model. In such a model 
one can accomodate many resonances coupled to interacting pairs of mesons. However, one should 
stress that the corresponding decay amplitudes are not unitary: unitarity is not preserved in the three- 
body decay channels and it is also violated in the two-body subchannels. Within the isobar model it 
is particularly difficult to distinguish the S -wave amplitudes from the non-resonant background terms. 
Their interference is often very strong which means that some two-body branching fractions, extracted 
from the data, could be unreliable. The isobar model is flexible but it has many free parameters (at least 
two fitted parameters for each amplitude component). For example, the Belle Collaboration has used 
40 fitted parameters in ref. |2| and the Babar Collaboration 43 free parameters in ref. Q]. 
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The construction of unitary three-body strong interaction amplitudes in a wide range of meson- 
meson effective masses is difficult. Therefore, as a first step, we incorporate two-body unitarity in our 
model of the D-decay amplitudes with final state interactions in the K^n* S- and f-waves and in the 
n + 7T~ S -wave. The branching fraction corresponding to a sum of these amplitudes is close to 80 % of 
the total branching fraction of the D — > K®n + n~ decay. 

2 Decay amplitudes 

The D° decays into K®n + n~ are analysed in the framework of the quasi-two body factorization ap- 
proach. The decay amplitude consists of the following 27 non-zero parts: seven allowed tree ampli- 
tudes, six doubly Cabibbo suppressed tree amplitudes, fourteen annihilation (W-exchange) amplitudes 
(7 allowed and 7 Cabibbo suppressed). The allowed amplitudes are generated by the quark c — > sud 
transitions and the doubly Cabibbo suppressed amplitudes correspond to c — > dus decays. Seven 
partial wave amplitudes are considered. We count S , P and D waves in the K^n* and the 7T + n~ subsys- 
tems and separately the P-wave amplitude for the G- parity violating <d — > n + n~ transition. There are 
numerous resonances contributing to different partial wave amplitudes. For example, in the K^n~ S- 
wave subchannel one can list 7T ( ,(800)~ or k~ and ^(1430)" resonances, in the P-wave one /f*(892)~, 
^1(1410) and A""(1680)~ and in the D-wave the #2(1430)" resonance. The same resonances but with 
the opposite charge are active in the K^n + subchannels. In the n + 7i~ subchannels one can enumerate 
the / (500) or <x, / (980) and / (1400) scalars, the p(770), p(1450) and w(782) vector states and the 
/2(1270) tensor resonance. 

The main part of the weak decay effective Hamiltonian is proportional to the following operator: 

O = %V' cs V ud hh, (1) 
V2 

where Gf is the Fermi coupling constant, V cs and V„d are the Cabibbo-Kobayashi-Maskawa matrix 
elements, the two quark currents are defined as j\ = (sc) V -a and j 2 = (ud) V -A, the symbol V - A 
means that the axial current is subtracted from the vector current. In the quasi-two body factorization 
approximation one can write: 

(K° s n-n + \ju 2 \D ) - (^tT^IDVI^IO) + (tTttI/^X^IO) + <0|/JD°><^7r7r + |./ 2 |0>, (2) 

where |0) denotes the vacuum state and j' x = (uc)v-a, f 2 - (sd)v-A are the quark currents obtained 
from the j\ , j 2 currents by the Fierz transformation. This form of the transition matrix element enables 
one to introduce the pion f„, kaon f K and D° f D decay constants. The corresponding expressions are: 

<;r + |j 2 |0> = if„ Px , (K° s \f 2 \0) = if K p K , (0\j[\D°) = -if D p D , (3) 

where p„, p K and p D are the pion, kaon and D meson four-momenta, respectively. 

2.1 Transition matrix elements 

The expression of the transition matrix elements can be simplified if we assume that the two strongly 
interacting mesons h 2 and hj, of momenta p 2 and pi form a resonance R in the final state. Then one 
can write: 

(h 2 ( P2 )h 3 ( P3 )\j \D°(p D )} = G Rhlh3 (s 23 )(R( P2 + p 3 )\j \D°(p D )), (4) 

where GRh 2 h 3 (s 23 ) with 523 = {pi + Pi) 1 is the vertex function describing the R decay into h 2 and h 3 
mesons. For example, let us consider the reaction D°(p D ) — > n + (pi)K^(p 2 )n^(p 3 ). If in the intermedi- 
ate state the #*(892) _ resonance is formed and decays into a K q tt~ pair then the transition form factor 
A D K " (ml) appears in the matrix element 

{R(P2 + P3)\j \D\pd)) = -2im K , e —^- p x A° K '~(ml) + 3 other terms. (5) 
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Here j - (sc)v-a, £ is the K*(892) polarization vector, m^. is the K* mass and po = P1+P2 + P3- The 
"3 other terms" do not contribute to the transition amplitude in eq. (|2). The vertex function G K ,-^o„- 
can be expressed in terms of the kaon-pion transition vector form factor Ff n (S23) 

G K '-K°A*2i) = e-(p 2 - p3)—^Ff°"'(s 23 ). (6) 

m K ,f K , 

In the above equation fx> is the K* decay constant. Assuming isospin symmetry one can equate the 
transition vector form factor Ff n (523) to the charged kaon to charged pion transition vector form 
factor Ff ^ (^23) calculated in ref. 0. 

The third term of eq. can also be simplified if the two hadrons, for example I12 and interact 
via a resonant state R. Then, similarly to eq. © we write 

(hi(pi)h 2 (p2)h 3 (p 3 )\f\0) = G Rh , h (s 23 )(h l (p 1 )R(p 2 + p 3 )l/|0>. (7) 

If, for example, h\ = K°, R = fo — > tt + tt~ and / = (M)v-a then the matrix element reads 

(K ( Pl )Mp2 + p 3 )|/|0) = -i " 1 * , ^ p D Ff fo (m 2 D ) + second term, (8) 

Pd 

k° f 9 

where F Q (m^) is the kaon to fo scalar transition form factor. The vertex function for the fo decay 
into jt + ti~ can be parametrized as 

G fo ^(s 23 )=X2Fo T+ *-(s 23 ), (9) 

where xi is a constant and Fi T+7I ~(s2 3 ) is the pion scalar form factor. Its functional form is taken from 
our B ± — > n ± n' 1 'n ± decay study |6). It preserves unitarity and groups together three scalar-isoscalar 
resonances / (500), / (980) and / (1400). 



2.2 Examples of allowed transition amplitudes with A'"/r final state interactions 



Using the assumptions introduced in subsection l2.1l one can derive formulae for the decay amplitudes 
in which the K®tt~ final state interactions are explicitely included. The S -wave amplitude reads: 

As = -^A iai Um 2 D - ml) F° K <> (m 2 „) F§°^(m 2 ). (10) 

Here, A\ = V* s V„d, a\ is the effective Wilson coefficient, F? (to 2 ) is the scalar D° to K*f transition 
form factor and F^"" (to 2 ) is the scalar K° to n~ transition form factor which depends on the K^n~ 
effective mass squared to 2 . The latter form factor can be taken equal to the K to n + transition scalar 
form factor F^ n (m 2 ) which has been evaluated in the study of B — > Kn + n~ decays Q. In this form 
factor both A"q(800) and ^(1430) resonances are included in a unitary way. The expression for the 
P-wave amplitude is: 

G F fn 1 1 (»in - to 2 )(to 2 - TO 2 ) „„»_ , «, , 

Ap = —f Aifli f^[TO 2 - to 2 + K —2 ^ *-] Aq (ml) Ff " (to 2 ), (11) 

^ JK" mi 

where to 2 is the K®n + effective mass squared and to 2 , denotes the n^if effective mass squared. The 
D-wave decay amplitude depends on the mass ntjc and the width Fk> of the resonance ^(1430): 



~A iai f n F DK i(m 2 _) 2 s , — . (12) 

2 m L K , - mi- im K ^F K ^ 



Here, F DK "-{m 2 _) is a combination of three types of the D to ,^(1430) transition form factors, g K -^ K 

is the decay coupling constant and the D(m 2 , to 2 ) is the D-wave angular distribution function. 

Derivation of other decay amplitudes proceeds in a quite similar way as shown above for the As , 
Ap and Ap amplitudes. One can notice that in our amplitude model the meson transition form factors 
play a very important role. 
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Fig. 1. Comparison of the Dalitz-plot projections for the Belle data [2] with the present model. 



3 Results 

The theoretical model outlined in the previous section has 28 free real parameters, most of them are 
unknown complex values of the meson-meson transition form factors appearing in the fourteen W- 
exchange amplitudes. The scalar and vector kaon-pion and scalar pion form factors are constrained 
using unitarity, analyticity and chiral symmetry. Fig. 1 shows that the model preliminary results are in 
fair agreement with the Belle Collaboration data. Also the total braching fraction is well reproduced 
0. 
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